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ABSTRACT 

We  describe  the  design,  implementation,  and  deployment  of  a  laser  slope  gauge  developed  at  the  Woods  Hole 
Oceanographic  Institution  for  the  purpose  of  studying  the  propagation  characteristics  of  capillary  ripples,  and  how  currents 
and  natural  slicks  on  the  ocean  surface  modify  ripple  spectra.  The  laser  slope  gauge  constitutes  a  nondisniptive  optical 
technique  for  determining  the  slope  spectrum  for  a  range  of  waves  with  wavelengths  between  2  mm  and  20  cm  using  both 
spatial  and  temporal  information.  Operation  of  the  sensor  and  data  acquisition  system  is  discussed  and  a  sample  data  record 
collected  in  the  Gulf  Stream  off  Cape  Hatteias,  NC  is  interpreted  and  analyzed. 


1.  INTRODUCTION 

The  study  of  short  wave  propagation  on  the  sea  surface!^  is  important  for  the  remote  sensing  of  subsurface  currents 
and  internal  waves  by  identifying  correlations  between  free  surface  hydrodynamics  and  features  observed  in  radar 
imageryf^*^!.  Multiple  sensor  measurements  gathered  during  the  Office  of  Naval  Research-sponsored  HIRES  I  experiment 
(the  pilot  experiment  of  the  high  resolution  advanced  research  initiative)  have  led  to  a  better  understanding  of  how  the 
collective  effects  of  wind  stress,  current  fields  and  surface  slicks  interact  to  produce  spatially  inhomogeneous  ripple  fields  on 
the  sea  surface. 

Central  to  the  field  experiment  was  the  deployment  of  a  scanning  laser  slope  gauge  (SLSG)  that  elicited  information 
about  the  propagation  of  ripples  for  studying  air/sea  interactions  and  how  slicks  on  the  ocean  surface  modify  the  wave 
spectrum,  e^recially  at  high  wavenumbcrsl^l.  The  instrument  was  designed  to  be  integrated  with  a  research  catamaran!^]  to 
acquire  unobstructed  measurements  of  wave  spectra  and/or  surface  roughness  as  the  platform  is  towed  at  a  rate  of  4  knots. 
The  dynamic  range  of  the  SLSG  is  consistent  with  gathering  spectral  information  of  waves  having  wavelengths  between  2 
mm  and  20  cm,  with  frequencies  up  to  64  Hz.  In  the  mode  of  operation  reported  herein,  only  frequency  spectra  are 
analyzed. 

The  catamaran  (referred  to  as  LADAS,  for  Laser  slope  gauge.  Acoustic  Dt^pler  current  profiler.  Acoustic  anemometer, 
Surface  sampler)  is  equipped  with  several  sensors.  Only  partial  results  obtained  from  the  SLSG  and  an  X-band  Doppler 
radar  ate  rqxxted  here. 


2.  SENSOR  DESIGN  &  DATA  COLLECTION  STRATEGY 


The  SLSG  consists  of  two  subassemblies,  namely,  a  laser  pod  and  an  optical  receiver.  Both  units  were  designed  to  be 
rigidly  attached  to  a  catamaran  with  the  laser  pod  deployed  below  the  waterline  and  the  receiver  above  the  waterline.  The 
basic  configuration  is  depicted  in  Figure  1.  As  the  catamaran  is  towed  through  the  water,  the  sensor  unobstnictively 
measures  the  local  slope  of  the  free  surface.  In  its  simplest  mode  of  operation,  the  laser  pod  transmits  a  stationary  vertical 
beam  that  penetrates  the  sea  surface.  In  this  mode,  the  apparent  frequency  spectrum  of  surface  slope  (up  to  83  Hz)  can  be 
measured.  In  its  more  complex  mode  the  laser  beam  is  scanned  to  extract  spatio-temporal  data  records  for  the  measurement 
of  full  3  dimensional  spectra  (with  spatial  and  temporal  resolution  of  10.5  cm'*  and  395  rad  s'*,  respectively).  We  will 
now  discuss  the  operating  principles  of  the  SLSG  and  show  results  obtained  in  the  apparent  frequency  mode. 
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Figure  lb.  Sensor  configuration  -  side  view 


As  illustrated  in  Figure  2,  the  beam  is  deflected  off  the  optical  axis  by  an  angular  amount  proportional  to  the  local 
slope  of  the  free  surface.  The  refracted  beam  is  focused  by  an  aspheric  lens  in  the  detector  head  so  that  its  image  at  the  focal 
plane  is  displaced  proportional  to  that  local  slope.  The  location  of  this  displaced  beam  is  detected  in  real-time  by  an 
orthogonal  set  of  four  photodiodes  in  the  detector.  The  centroiding  algorithm  relies  on  the  individual  measurements  from 
each  of  the  photodiodes  to  infer  the  angular  deflection  of  the  beam,  and  thus,  the  x  and  y  components  of  surface  slope.  As 
the  catamaran  is  towed,  the  collected  time-series  represents  the  variation  of  ripple  slope  with  distance  -  which  is  the  spatial 
derivative  of  ripple  amplitude.  Inertial  sensors  onboard  the  catamaran  allow  for  the  subtraction  of  the  platform  velocity  to 
compute  the  wavelength  of  the  ripples. 


sin'*(1.33sina)-a  =  (p  ==>  (p«=0.33a 


Figure  2.  Beam  geometry  at  the  air/sea  interface 


2.1  The  Laser  Pod 

The  laser  pod,  shown  in  Figure  3,  contains  a  10  milliwatt  cw  helium-neon  laser,  a  two-axis  galvo  scanner,  two  fold 
mirrors,  and  a  custom  aspheric  lens.  All  the  components,  except  for  the  lens,  are  secured  to  a  common  tray  which  is 
attached  to  one  of  the  end  caps.  The  tray  facilitates  the  convenience  of  performing  the  optical  alignment  outside  of  the 
tube,  in  addition  to  keeping  the  alignment  intact  whenever  the  tray  is  removed  from  the  tube  for  diagnostic  purposes. 

The  He-Ne  laser  emits  a  0.8  mm  diameter  TEMqo  beam  measured  at  the  1/e^  irradiance  points.  The  two-axis  scanner 
employs  two  limited  rotation  servo  motors  (galvonometers)  configured  in  an  orthogonal  set.  Each  galvo  controls  the 
movement  of  a  small  minor  to  perform  beam  steering  in  the  x-  and  y-directions.  A  capacitive  rotation  sensor  in  each 
device  provides  closed  loop  control  over  the  angular  position  of  the  minors. 

Beam  {Hopagation  through  the  laser  pod  can  be  described  by  the  following  sequence  of  events.  Fold  mirror  Ml 
redirects  the  light  emerging  from  the  laser  by  90°  up  to  the  scanner.  At  the  scanner,  the  beam  strikes  the  x-galvo  minor 
first  and  bounces  a  few  millimeters  away  to  the  larger  y-galvo  mirror.  After  the  y-galvo  minor,  the  beam  has  been  slewed 
in  two  dimensions.  The  resulting  scan  pattern  is  propagated  down  the  axis  of  the  laser  pod  to  a  large  fold  mirror  M2  which 
redirects  the  light  90°  up  to  the  asphere  lens  LI.  The  lens  collimates  the  scan  pattern  so  that  all  off-axis  cays  become 


parallel  once  they  emerge  from  the  lens.  The  lens  forms  the  exit  aperture  of  the  laser  pod  and  was  designed  to  operate  with 
air  on  one  side  and  sea  water  on  the  other. 


The  laser  pod's  optical  axis  is  co-aligned  with  the  detector  head  so  that  any  angular  offset  of  the  beam  from  a  parallel 
ray  is  solely  an  artifact  of  the  intervening  free  surface  interface.  After  penetrating  the  free  surface,  the  beam  is  captured  by 
the  receiver  where  its  angular  deflection  is  measured  in  terms  of  the  spot's  centroid  on  the  horizontal  x-y  plane. 


Section  A-A* 
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Figure  3.  Opto-mechanical  layout  of  the  laser  pod 


2.2  The  Detector  Head 

A  schematic  of  the  detector  head  is  shown  in  Figure  4.  Along  the  forward  direction  of  the  optical  path,  components  in 
the  head  unit  include  a  bandpass  filter  (A),  a  plano-convex  lens  (B),  a  diffuser  plate  (Q,  an  image  plane  reflector  plate  (D), 
and  a  printed  circuit  board  (E)  containing  four  silicon  photo-diodes  with  10  nm  interference  filters  and  analog  signal 
conditioning  electronics. 


The  interference  filters  are  centered  about  the  633  nm  He-Ne  line  for  strong  ambient  light  rejection.  The  diffuser  plate, 
positioned  at  the  focal  plane  of  the  lens  L2,  acts  as  a  quasi-Lambcrtian  surface  operating  in  transmission.  Its  purpose  is  to 
forward  scatter  the  laser  radiation  into  a  broad  Gaussian  profile.  Quad  cell  detection  at  a  distance  of  12  cm  above  the  diffuser 
is  used  to  cenuoid  the  location  of  the  beam.  The  centroiding  algorithm  yields  measurements  of  the  x-  and  y-componcnts  of 
surface  slope. 


Figure  4.  Schematic  of  the  Detector  Head 

23  Sensor  Control,  Data  Acquisition  &  Modes  of  Operation 

A  functional  block  diagram  showing  the  electrical  interfaces  between  the  sensor,  control  electronics,  and  data 
acquisition  system  is  illustrated  in  Figure  S.  In  addition  to  the  laser  pod  and  detector  head,  the  overall  system  hardware 
includes: 

•  a  Goieral  Scanning  DE'2(XX)™  scan  controller  which  drives  the  galvos  in  the  laser  pod,  controlled  by 
serial  line  (E): 

•  two  programmable  microcontrollers  [WHOI-SOOO's]  that  read  n  data  points  in  response  to  a  TTL  trigger 
supplied  by  the  DE-2000  via  data  line  {D} ,  (the  number  n  is  a  function  of  the  sampling  scheme); 

•  a  25  MHz  386  computer  with  two  Metrabyte  PlO-12  parallel  I/O  cards,  a  3COM  Ethernet  adapter, 
and  serial  port; 

•  a  DEC  5000  computer  on  the  Ethernet  line  (G) ; 

•  a  Compaq  portable  on  the  serial  lines  (E)  and  (F). 


WHOI 3000 
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Figure  5.  Functional  block  diagram  of  the  sensor  system 


The  386  receives  messages  at  the  serial  line  to  initiate  sampling.  Data  collection  activity  proceeds  until  another 
message  requests  the  386  to  stop  sampling.  The  DE-2000  is  programmed  to  execute  scan  sequences  for  spatio-temporal 
sampling  and  non-scanning  modes.  During  the  square  pauem  sampling  mode,  the  scanner  writes  out  a  square  with  10  cm 
sides.  At  the  start  of  each  square,  the  DE-2000  emits  a  TTL  pulse.  The  WHOI-3000's  then  begin  A/D  conversions  and 
each  WHOI-3000  samples  a  light  intensity  value  corresponding  to  its  reqrective  x-  or  y-channel.  Each  will  accumulate  129 
pairs  of  readings  in  about  10  ms.  This  corresponds  to  32  pairs  per  side  of  the  square.  Each  pair  represents  the  light 
intensities  at  the  two  photodiodes  for  one  component  of  slope. 


3.  DATA  ANALYSIS 

In  this  Section,  a  sample  time  series  of  2-D  surface  slope  acquired  by  the  laser  slope  gauge  is  analyzed.  The  statistical 
and  spectral  properties  of  the  surface  slope  data  are  related  to  sensor  parameters  and  the  efficacy  of  the  sensor  to  measure 
ocean  surface  toughness  is  discussed.  First,  we  begin  with  a  description  of  the  pre-processing  p^ormed  on  the  stored  data 
to  prepare  it  for  data  analysis. 


3.1  Pre-processing  of  the  Raw  Data 


The  spot  position  on  the  image  plane  as  a  function  of  time  was  constructed  using  the  four-channel  time-series  of 
intensities.  The  (x,y)  coordinates  of  the  spot  position  were  calculated  using  a  component-pair  centroiding  algorithm: 


fore  /  aft  coordinate 
transverse  coordinate 


which  produces  a  normalized  position  in  the  local  coordinates  of  the  image  plane  with  limits  -1  S  x  ^  1  and  -1  ^  y  1,  as 
illustrated  in  Hgure  7.  When  the  intensities  at  the  four  detectors  are  equal,  the  computed  cenUDid  of  the  spot  is  (0,0),  i.e. 
the  spot  is  aligned  with  the  optical  axis,  corresponding  to  a  zero  surface  slope  condition.  Given  the  geometry  of  the 


refracted  beam  (see  Figure  2)  and  the  optical  ray  trace  of  the  head  unit  (see  Figure  4),  the  position  of  the  spot  on  the  image 
plane  is  direcdy  related  to  the  local  inclination  of  the  sea  surface.  For  small  angles  (^  13  degrees)  of  beam  deflection,  the 
spot  position  on  the  image  plane  is  related  to  the  surface  slope  in  the  x-  and  y-dircctions  by  the  relationship: 

ax  *=  3  tan’^  (Cx)  ==>  fore/aft  surface  slope 

ay  -  3  tan'*  (Cy)  ==>  transverse  surface  slope 

where  C  is  a  constant  defined  by  the  ratio  of  the  spatial  separation  of  the  detectors,  C.  from  the  optical  axis  and  the  focal 
length,/,  of  the  lens  in  the  receiver.  For  our  system,  ^  =  8.85  cm  and /  =  43.18  cm  which  makes  C  =  0.205.  Hence,  the 
laser  slope  gauge  is  capable  of  measuring  surface  slopes  of  ±  34.8*. 

Upon  completion  of  the  pre-processing,  three  time-series  are  generated  from  the  raw  data: 

•  Four  channels  of  detector  intensities  Ia.  Ib.  IC  and  *01 

•  Centroided  spot  position  (x,y)  on  the  image  plane; 

•  2-D  surface  slope  ax  and  ay. 

The  detector  intensities  take  on  values  between  0  and  1023  as  a  result  of  our  10-bit  AJD  conversion,  and  the  normalized 
spot  position  and  surface  slope  have  the  limits  -1  ^  x,y  5  -tl  and  -35*  ^  ax,y  5  +35*,  respectively. 
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Figure  7.  Normalized  spot  position  on  the  image  plane 


3.2  Description  of  the  HIRES  I  Data  Record 

During  the  HIRES  I  Held  experiment  off  the  coast  of  Cape  Hatteras,  the  laser  slope  gauge  was  operated  in  the  non¬ 
scanning  mode  for  the  measurement  of  ocean  surface  roughness.  The  data  record  reported  here  was  collected  on  the  16th  of 
September  beginning  at  10:09  AM  local  time  and  lasts  for  approximately  10  minutes.  According  to  the  ship's  log  and  data 
acquired  by  other  sensors,  the  laser  slope  gauge  data  should  contain  features  associated  with  a  current-field  boundary.  Based 
on  satellite  thermal  imagery,  it  is  believed  that  this  current  field  may  be  associated  with  an  eddy  along  the  wekem  edge  of 
the  Gulf  Stream. 


The  catamaran  was  being  lowed  at  4  knots  («=  2  m/s)  and  data  was  sampled  at  167  Hz.  Given  this  tow  velocity  and 
sampling  rate,  the  platform  was  translated  approximately  1,2  cm  between  samples.  The  spatial  resolution  of  the  sensor  is 
governed  by  the  1  mm  diameter  of  the  laser  beam. 

It  is  important  to  note  here  that  the  measurements  are  somewhat  difficult  to  interpret  owing  to  their  mixed  temporal 
and  spatial  nature.  For  a  fixed  laser  slope  gauge,  time-series  measurements  can  be  analyzed  by  Fast  Fourier  Transforms 
(FFT)  to  produce  frequency  spectra.  For  a  scanning  laser  slope  gauge,  FFT  analysis  directly  yields  wavenumber  spectra. 
The  measurements  obtained  in  the  September  field  experiment  are  a  hybrid  of  these  techniques  due  to  the  fact  that  the  fixed 
laser  was  translated  at  about  2  m/s.  This  rate  is  large  enough  so  that  some  waves  (those  with  phase  velocities  less  than  2 
m/s)  were  overtaken  by  the  gauge,  and  others  were  not  This  has  the  effect  of  aliasing  waves  into  improper  frequency  (or 
wavenumber)  bins.  To  avoid  errors  resulting  from  this  effect,  the  analysis  described  below  should  not  be  taken  to  literally 
represent  either  frequency  or  wavenumber  spectra.  It  is  intended  to  represent  surface  roughness  via  the  frequency 
representation  of  slope  variance. 

3.3  Statistical  Analysis 

Evaluating  the  statistical  properties  of  the  sample  data  record  reveals  important  information  about  the  performance  of 
the  sensor.  Parameters  such  as  sensor  bias,  dynamic  range,  and  scale  factor,  for  example,  can  be  estimated  from  the 
computed  mean,  standard  deviation,  and  covariance  of  the  detector  intensity  and  spot  position  time-series. 

If  the  laser  pod  and  detector  head  were  perfectly  aligned  and  the  responsivity  of  the  detectors  were  equal,  we  would 
expect  identical  values  of  mean  detector  intensity  over  the  10  minute  data  record  for  the  four  optical  channels,  i.e. 


hi  practice,  however,  this  is  rarely  the  case.  At  the  time  this  data  was  collected,  the  mean  values  of  detector  intensity  were 

1^=430  r^=470 

10  =  449  r0,=  431 

which  yields  a  mean  centroided  spot  position  in  the  normalized  coordinate  frame  of 

x  = -4.165  X  10*2 
y=  6.067  X  10*2 

corresponding  to  a  surface  slope  bias  of 

a^  =  -1.47' 

ay  =  2.14‘ 

This  surface  slope  bias  is  related  to  the  cumulative  errors  of  both  nonuniform  responsivity  and  optical  misalignment, 
resulting  in  an  effective  off  set  from  the  optical  axis  given  by 

r=  V(i)2+  (7)2 

=  .074 

where  r  is  the  radial  distance  away  from  the  optical  axis,  in  the  normalized  coordinate  frame,  whose  real  distance  is 

R  =  rC 


=  .65  cm 


The  effecUvc  origin  of  the  image  plane  is  therefore  6.5  mm  aft  and  to  port  of  the  real  optical  axis.  This  brings 
asymmetry  to  the  dynamic  range  of  our  sensor.  Rather  than  having  the  design  measurement  limits  of  ±34.8*  in  x  and  y 
surface  slope,  we  instead  have  [-33.3*,  +36.2*]^  and  [-36.9*,  +32.6*]y.  An  adjustment  to  our  scale  factor  can  be  made  that 
accounts  for  the  bias. 

Ox  =  1 .47  +  (3  tan'  ^  (.205x)]  ==>  fore  /  aft  surface  slope 

ay  =  -2.14  +  [3  tan'^  (.205y)]  =>  transverse  surface  slope 

The  covariance  of  the  spot  coordinates  a\y  also  serves  as  a  useful  diagnostic.  Figure  8  is  a  scatter  plot  in  ax  :  Oy 
space  with  a  2a  covariance  ellipsoid  superposed  on  the  data.  Each  point  on  this  plot  corresponds  to  the  standard  deviation 
in  the  x-  and  y-spot  coordinates  simultaneously  computed  over  I  second  intervals  containing  167  samples. 
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Figure  8.  Covariance  of  the  x,y  spot  coordinates  on  the  image  plane 


Since  the  standard  deviation  of  spot  position  is  related  to  the  magnitude  of  surface  roughness,  the  +45*  orientation  of 
the  ellipsoid  is  expected  for  these  positively  correlated  parameters.  In  other  words,  when  high  surface  roughness  is 
encountned  along  the  x-axis,  high  surface  roughness  should  (on  average)  be  encountered  along  the  y-axis  as  well.  In  the 
limit  of  a  perfectly  correlated  distribution,  the  ellipsoid  would  collapse  to  a  line.  Due  to  the  nature  of  the  process  and 
measurement  scheme  we  do  not  anticipate  a  perfectly  correlated  distribution  in  Cx :  ay  space.  What  we  do  expect,  however, 
is  a  smooth  distribution  of  the  covariance  without  abrupt  transitions  or  evidence  of  multimodal  behavior.  The  laser  slope 
gauge  does  indeed  exhibit  this  smoothness  as  witnessed  in  Figure  8. 


As  a  first-cut  procedure  for  scanning  the  data  for  surface  roughness  features,  it  is  convenient  to  plot  Cx(0  and  cry(t). 
Figure  9  shows  such  a  plot  for  the  time  interval  of  400  <  t  5  500  seconds.  The  standard  deviations  Ox  and  Oy  were 
computed  in  the  same  manner  as  the  results  of  Figure  8,  i.e.  at  1  second  intervals  using  167  samples. 
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Figure  9.  Time*history  of  the  x-  and  y>spot  coordinate  standard  deviations 

Figure  9  contains  three  distinct  features,  namely,  a  wide  lobe  of  high  Ox  and  Oy  centered  near  t  =  430  sec,  a  region  of 
relatively  low  Ox  and  Cy  lasting  for  15  seconds  around  t  =  450  sec,  followed  by  a  strong  and  highly  correlated  modulation 
beginning  at  t  =  460  sec  and  lasting  for  35  seconds.  Based  on  our  tow  velocity  of  4  knots,  these  three  features  have  a  spatial 
extent  [in  the  direction  of  platform  heading]  of  1.6  meters,  2.4  meters,  and  7.0  meters,  respectively.  The  period  of  the 
modulating  feature  starts  out  at  approximately  8  seconds  and  then  "bunches"  towards  a  2  second  period  before  diminishing. 
The  physical  interpretation  of  this  feature  may  be  related  to  bands  of  roughness,  which  will  be  explored  further  in  Section 
3.4. 


3.4  Spectral  Analysis 

In  this  section,  short-time  Fourier  transform  (STFT)  representations  of  the  sample  data  record  arc  presented  to 
characterize  the  change  in  frequency  content  of  the  surface  roughness  over  time  as  the  catamaran  is  towed  through  the  water. 
Attention  will  be  focused  on  the  region  of  data  suspected  of  containing  peculiar  surface  features  found  by  statistical 
methods.  The  sensor’s  ability  to  detect  and  quantify  roughness  variations  on  the  sea  surface  are  addressed  by  reviewing 
video  and  radar  backscattcr  data  collected  over  the  same  time  domain. 


The  underlying  idea  for  the  STFT  is  lo  take  a  discrete  Fourier  transform  (DFT)  in  the  neighborhood  of  each  time 
sample  of  our  sequence.  The  resulting  STFT  representation  will  be  related  to  both  the  time-  and  frequency-domain 
properties  of  the  observed  surface  roughness.  The  signal  processing  techniques  employed  are  an  extension  of  the  basic 
Fourier  transform  concepts  for  a  discrete-time  sequence  x(n).  The  discrete  STFT  of  x(n)  is  defined  by  the  expression^'®) 

m  =  oo 

X(njc)  =  x(m)  w(n  -  m) 

m  = 

where  n  and  k  are  discrete  variables  denoting  time  and  frequency,  respectively,  w(n  -  m)  is  the  analysis  window,  and  N  is  the 
frequency  sampling  factor. 

Figure  10  gives  the  STFT  result  obtained  with  a  64-point  Hamming  window  with  double  zero  padding.  The  plotted 
result  has  a  linear  frequency  scale  of  0  -  83  Hz  (increasing  frequeiKy  to  the  left),  and  is  associated  with  the  time  interval  380 
^  t  S  440  seconds  (increasing  lime  is  from  front  to  back).  Thresholding  was  applied  for  suppression  of  background  noise. 
No  additional  information  of  substance  was  revealed  from  the  use  of  a  128-point  window  excqit  the  realization  that  the  data 
is  very  wide  band. 


Figure  10.  STFT  result  using  a  64-point  Hamming  window 


Figure  11  Comparison  of  STFT  with  spot  coordinate  standard  deviation 

and  X-band  Doppler  radar  time-series 


Figure  1 1  shows  a  comparison  of  the  STFT  {contour  plot)  with  the  time-history  of  standard  deviation  and  X-band 
Doppler  radar  return  power  over  the  time  frame  380  <  t  ^  560  seconds.  The  X-band  Doppler  measurements  were  made 
independently  from  a  mast-mounted  Doppler  scatterometer  operated  by  Dr.  David  Lyzenga  of  ERIM.  X-band  Doppler  return 
power  is  proportional  to  surface  roughness  and  thus  serves  as  an  independent  measurement  of  roughness.  Perfect 
agreement,  however,  is  not  expected  because  of  two  reasons: 

1)  backscatter  power  is  a  strong  function  of  angle  of  incidence  and  the  catamaran  is  constantly  undergoing  pitching 
motions  (because  the  scatterometer  was  not  mounted  on  a  stabilized  platform,  angular  fluctuations  due  to  pitch 
directly  affect  angle  of  incidence),  and 

2)  the  measurement  regions  of  the  radar  and  laser  slope  gauge  are  physically  separated  from  each  other  (by  5  m 
spatially,  and  2.5  seconds  temporally)  and  therefore  some  surface  features  detected  by  one  sensor  are  not 
necessarily  detected  by  the  other. 

A  reasonable  correlation  exists  between  the  STFT  and  the  X-band  radar  data.  This  correlation  suggests  the  existence  of 
some  features  that  appear  on  the  surface  long  enough  so  as  to  be  recorded  by  both  sensors.  These  features  might  include 
wind-stress  induced  surface  roughness  or  current  boundary  surface  expressions.  A  strong  correlation  exists  between  features 
observed  in  the  STFT  and  standard  deviation  plots.  This  is  to  be  expected  since  the  STFT  surface,  if  integrated  along  the 
frequency  axis,  yields  the  mean-square-slope,  a  quantity  akin  to  standard  deviation  as  defmed  and  computed  above. 

The  correlation  between  the  STFT,  standard  deviation,  and  video  signal  is  less  quantifiable  since  the  video  is 
unprocessed.  All  that  can  be  assessed  is  the  presence  (or  absence)  of  specular  reflections  from  facets  on  the  sea  surface. 
Modulations,  or  occunences  of  alternately  rough  and  smooth  surfaces,  were  evident  throughout  the  video  for  the  times 
analyzed  above. 

4.  CONCLUSION 

The  Laser  Slope  Gauge  was  successful  in  measuring  surface  roughness  over  a  variety  of  sea-state  conditions  during  the 
HIRES  I  field  experiment  off  the  coast  of  North  Carolina.  The  sensor  was  operated  exclusively  in  a  non-scanning  mode  for 
the  purpose  of  understanding  the  modification  of  surface  roughness  in  the  presence  of  current  Held  boundaries  and  surface 
water  temperature  fronts.  The  utility  of  the  LADAS  catamaran  to  perform  "ground  truthing"  for  airborne  radar  platforms 
was  also  demonstrated  during  the  experiment 

Data  sets  collected  by  the  Laser  Slope  Gauge  have  been  processed  and  analyzed  in  conjunction  with  data  acquired  by  the 
other  sensors  on  the  LADAS  catamaran.  The  processing  procedures  included:  1)  centroiding  the  spot  position  on  the  image 
plane  and  2)  statistical  and  spectral  analysis  to  determine  instrument  bias  and  dynamic  range.  The  collected  data  has  been 
compared  in  both  a  single-dimensional  format  (e.g.  standard  deviation  as  a  function  of  time)  and  in  a  multi-dimensional 
representation  (e.g.  STFT ,  in  both  time  and  frequency)  to  independent  measurements  of  surface  roughness. 

The  sensor's  scanning  mode  of  operation  has  been  tested  in  the  laboratory  and  used  for  the  first  time  in  the  field  during 
an  experiment  in  May  of  1992  off  the  Massachusetts  coast  in  Martha's  Vineyard  Sound.  One  of  the  major  thrusts  of  this 
expCTiment  was  to  gather  directional  spectra  [i.e.  S(kx.  ky,  ©)  ]  of  capillary-gravity  waves  in  the  presence  of  natural  slicks. 
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